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Abstract 

In this paper, we present an algorithm of simple exponential growth called 
COPOMATRIX for determining the copositivity of a real symmetric matrix. 
The core of this algorithm is a decomposition theorem, which is used to deal 
with simplicial subdivision of T~ = {y G A m \(3 T y < 0} on the standard 
simplex A m , where each component of the vector is -1, or 1. 
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1. Introduction 

QUESTION 1 Let A be a given n x n real symmetric matrix, M + be 
the set of nonnegative real numbers, and 

Q(X) = X T AX, X ^ 

be a quadratic form. What conditions should A satisfy for [VX G R" , Q(X) > 
0(>0)]? 

If [VX G R", Q(X) > (> 0)], then the quadratic form Q(X) is called a 
(strictly) copositive quadratic form and the corresponding matrix A is called 
a (strictly) copositive matrix. 
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Copositive matrices have numerous applications in diverse fields of ap- 
plied mathematics, especially in mathematical programming and graph the- 
ory (see @, S, 0, Ell, Q E H S H 0H). Therefore copositivity has 
been studied thoro^lily since 1950s (seeQJJ Q 0, 0, 0, B, [33 



In general, it is an NP-complete problem to determine whether a given 

This means that every 



n x n symmetric matrix is not copositive [37|, 138 



algorithm that solves the problem, in the worst case, will require at least 
an exponential number of operations, unless P=NP. For that reason, it is 
still valuable for the existence of so many incomplete algorithms discussing 
some special kinds of matrices (see @, |, 0, S, H El 0, H H ) • For 
small values of n(< 6), some necessary and sufficient conditions have been 
constructed (see [H, 14, 30, 0]). From another viewpoint, QUESTION 
1 is a typical real quantifier elimination problem 0, 0, 0, 35, 40, 44, 
45], which can be solved by standard tools of real quantifier elimination 
e.g., using CAD) Isl Eg . |47| . Thus, there is a complete algorithm for 

determining copositive matrices theoretically. Unfortunately, this algorithm 
is not efficient in practice for the CAD algorithm is of doubly exponential 
time complexity (see 0, 0, B])- In this paper, we will construct a complete 
algorithm with singly exponential time bound. 

The standard simplex A m (m > 2) is defined as the following set 



{(yi,- 



,Vm) 



2/1 H 1" Vm = 1, Vl > 0, ■ ■ ■ , Vm > 0}. 



It is well known that the dimension of A m is m 



e m , namely, e x = (1, 0, 



Let A G 



, '-"mi 
Draxn 



Of 



1. Denote the vertices of 
m = (0,0,...,lf. 



be symmetric and be partitioned as 



A 



\Oti 



an a 
a A 2 



> 



Define B = anA 2 — aa T . It is easy to see the following facts (cf. [l|) 

1. If aiij > 0, i = 2, . . . , n, then A is (strictly) copositive <^=^ an > 
0) and A 2 is (strictly) copositive. 

2. If at least one of an is negative, then we need only to focus on the 
set of points T~ = {y G A n ^i\ a T y < 0}. It is well known that T~ is a 
convex polytope on A n _! (see [l(]). The polytope T~ can be subdivided into 



2 



the simplices Si,...,S p , that is, 



T~ = {JSi, int(S' 4 )f|mt(S , ,) = 0,for^ 



i=i 



where int(Si) denotes the interior of simplex Si. The coordinates of the ver- 
tices that span the simplex Si constitute a matrix denoted as W{. Andersson 
et al. p. 23) proved the following results. 



Lemma 1.1. (a) A is copositive iff a u > and A 2 , W^BW U W?BW P 
are all copositive. 

(b) A is strictly copositive iff an > and A 2 , Wl BW%, . 
are all strictly copositive. 



W^BW P 



In order to formulate the algorithm of Lemma ll.l[ we first consider how to 



{y e A r 



a 



y<o}- 



obtain the simplicial subdivision of the polytope T" 
For small values of n(< 6), Andersson et al.[l] and Yang and Li 49] give the 
simplicial subdivision of T~ . However, they do not provide a procedure for a 
simplicial subdivision of T~ for arbitrary values of n. We propose a simplicial 
subdivision of T" for all values of n, and consequently construct a complete 
algorithm for determining the copositivity of an n x n matrix. 

We will adopt a flexible approach. Rather than subdivide T~ into sim- 
plices (of course our method is also valid for subdividing T~ into simplices), 
we first transform the matrix A into the following matrix called A. 

Let a = («i2, • • • , «i„) T and D = diag(c?2, • • • , d n ), where 



Then 



.4 



1 
D 



A 



f\aii\, 

1 
D 



if oc u = 0; 
if an 0. 



a DA 2 D 



(1) 



where a = (sign(ai 2 ), • • • , sign(«i n )) T . 

Obviously, A is (strictly) copositive 
Lemma [1.11 to A. Let 



'A is (strictly) copositive. Apply 



01 



sign(ai 2 ; 



sign(ai n j 
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Thus we just need to subdivide T into simplices, where 

f- = { y e A„_! . . . , &_!)!/ < o, A e {-1, o, 1}}. 

Next we make further simplification: separate -1,0,1 from . . . , /3 n _i, 
namely let 

An = • • • = = 1, A* = • • • = bt = -1, CI = ■ ■ ■ = f3 Cr = 0. 

{ai, . . . , a s , ci, . . . , c r } = {1, . . . , n - 1}, 

r, s,t>0, t > 1, r + s + t = n-l. 

In geometry it is easy to see that the convex polytope T~ is the convex hull 
of its surface S~ and its vertices e Cl , . . . , e Cr , that is, 

T~ = conv{e cl , . . . , e Cr , S 1- }. (2) 

s ~ = {(yu ■ ■ -,yn-i) T e A n _i|y 0l h hy^-y^ Vh < 0, ^ 

If the simplicial subdivision of S 1 " is known, the simplicial subdivision 
of T~ is directly obtained by d2J). So we just need to study the simplicial 
subdivision of the polytope S~ . 

2. A simplicial subdivision algorithm for the convex polytope S~ 

2.1. Fundamental notations 

The notation S~ is simple, but it can not reveal the information of convex 
polytopes. In order to simplify the descriptions, we will introduce a new 
notation, which is fundamental to our study. 

Definition 2.1. Suppose that two sequences of positive integers [a 1; a%, . . . , a s ], 
■ ■ -A] satisfy 

{ai,...,a s , &i, . . . , b t } C {1, 2, . . . , m}, s > 0, t > 1, m > s + 1 > 2, 

where all of s + t elements of {ax, . . . ,a s , b\, . . . , b t } are distinct. Then the 
notation [[ai, a 2 , . . . , a a ], [6i, b 2 , ■ ■ ■ , bt]} m is defined as the polytope S~ in (E]). 
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For example, let us compare the polytope [[2,3], [5]] 5 and the polytope 
[ [2 , 3] , [5] ] 6 . [ [2 , 3] , [5] ] 5 denotes the polytope 

{(yi, • • • , 2/s) T e A 5 \y 2 + y 3 - 2/5 < 0, (y 2 , y 3 , y->) T e A 3 }. 

Here (2/2, 2/3, ys) T £ A 3 implies that yi = 0, 2/4 — 0. [[2,3], [5]] 6 indicates the 
polytope 

{(2/1, ■ ■ • , 2/e) T e Aeli/a + 2/3 — 2/5 < 0, (3/2, 2/3, 2/s) T e A 3 }. 

Here (2/2, 2/3, 2/s) T £ A 3 implies that 2/1 = 0,2/4 = 0,2/6 = 0. It is clear that 
[[2, 3], [5]] 5 and [[2, 3], [5]] 6 are congruent, although they are sets in simplices 
of different dimensions. 

For < k < m — 1, the polytope is defined as 

L k = {{Vu Vrnf G A m |yi H h y k - Vk+i y m < 0}- 

is written as [[1, . . . , k], [k + 1, . . . , m]] m by the notation of Definition 
12.11 is a special case of S~ , but this notation is more convenient for our 
analysis. 

In the following we will study the basic geometric properties of convex 
polytope [[ai, a 2 , . . . , a s ], [h, b 2 , . . . , b t }} m . 

2.2. Geometric properties of S~ 

Let ei, . . . , e m be vertices of the standard simplex A m , and My = (e^ + 
e 3 -)/2 be the midpoint of the line segment e^e]. 

The following result is stated in [l[ without proof. For completeness, we 
give a proof. 

Lemma 2.1. [1] Given a convex polytope , then all of its vertices are 

V = {e k+1 , . . . , e m , Mjj, £ = 1, 2, . . . , k, j = k + 1, . . . , m}. 
The number of the vertices is |V| = (k + l)(m — fc). 

Proof. Note that the convex polytope is obtained by cutting the stan- 
dard simplex A m with the hyperplane 

L=o : Vi H h 2/fc - 2/fc+i 2/m = 0. 
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Therefore the vertices of the polytope come from two parts: one part is 
vertices of A m , that is, {e^+i, . . . , e m }; while the other part is the intersection 
points of the hyperplane L =0 and the edges of standard simplex A m . 

First consider the intersection point of L =0 and the edge aei+beu+i (a, b > 
0, a + b = 1). Substitute ae± + bek+i into the following equations, 



Therefore, the solutions are a = 1/2, b = 1/2, namely, the intersection point 
is Mi )fc+ i. 

In the same way, we get all intersection points of L = o and the edges of 
A m . They are {Mjj, i — 1, 2, . . . , k, j = k + 1, . . . , to.}. 

Hence the number of vertices of is \V\ = m — k + A;(m — = (k + 
1)(to — fc). 

Likewise, we can prove the following lemma. 

Lemma 2.2. Given a convex polytope [[ai, a 2 , . . . , a s ), [bi, b 2 , ■ ■ ■ , b t ]) m , then 
all of its vertices are 

V = {e bl , . . . ,e bt , M aitbj ,i = 1,2, ...,s, j = 1, . . . ,t}. 

The number of the vertices is |V| = (s + l)t. 

We see that the polytope [[ai, ... , a s ], [6i, 62, • • • , &t]]m an d the polytope 
L fc are similar in many respects, which will be further discussed. 

Lemma 2.3. The convex polytope is simplicial iff k — 0, or k — to — 1. 

Proof. When k = 0, L = A m is simplicial. 

When k = m — 1, consider the convex polytope 

L m-i '■= {{Vu ■ ■ ■ > l/m) T e A m | H h y m _i - y m < 0}. 



By Lemma [27TI we know that all vertices of L^ n _ 1 are {e m , Mi jm , i = 1, 2, . . . , 
to — 1}. Obviously all the vectors of {M it7n — e m , i = 1, 2, . . . , to — 1} are 
linearly independent, so is simplicial. 

Conversely, we know that the dimension of the polytope is to — 1. 



If k 7^ 0, to — 1, then by Lemma I2.1[ the number of the vertices of L fc is 
(k + l)(m — k) 7^ to, so is not simplicial. 



l/i H \-Vk- Vk+i Urn 



0, 2/1 H hi/, 



1 
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Lemma 2.4. The convex polytope [[a\, a 2 , . . . , a s ], [61, 62, ... , b t }}m (here the 
vertices are obtained by Lemma 12. 2 p is simplicial iff s = 0, or t = 1. 

Lemma 2.5. The dimension of the polytope [[a\, 02, ... , a s ], [61, 62, ... , bt]] m 
is (s + t- 1). 

If the polytope [[ai, 02, ... , a s ], [61, 6 2 , . . . , b t }]m is not a simplex, we will 
subdivide it into simplices. 

Lemma 2.6. If the polytope [[ai, a 2 , . . . , a s ], [61, 6 2 , ■ ■ ■ , b t )) m is not a sim- 
plex, then there are only two (s + t — 2)-dimensional surfaces that do not 
include the vertex M 0l) 6 x . They are 

[[a 2 , ■■■,a s ], [h, 6 2 , ••• , &t]]m, ■■■,a 8 ], [b 2 , • • • , &t]] m . 

(obtained by deleting a\, b\ from array [[ai, a 2 , . . . , a s ], [61, b 2 , . . . , b t ]}m respec- 
tively) 

Proof. All the (s + t — 2) -dimensional surfaces of the convex polytope 
[[ax, a 2 , . . . , a s ], [h, b 2 ,..., b t )] m are obviously 



[\ai,a 2 , . . 


■ ,a s ], [61,62, • 


• • , 6t]] m , 


[[oi,o 2 , • • 


■ ,a>s], [61,62,- 


• - , 6t]] m , 


. . . , 

[[a u a 3 , . . 


• ,a s ], [61,62, . 


• • ,b t )) m 



(where the notation [[ai, a 2 , ... , a s ], [61, 62, ... , 6j]] m is the polytope with ai 
removed) and 

{{yi, ■ ■ ■ ,y m ) T e A m |y ai h \-y as -yh 2/^ = 0, (y ai , . . . ,y as , 

Ubi, ■ ■ -,yb t ) T e A s+ t}. 

That makes s + t + 1 (s + t — 2)-dimensional surfaces in all. We can verify 
that only 

[[a 2 , . . . , a s ], [61, 6 2 , ... , b t ]] m , [[ai, ...,a s ], [b 2 , . . . , b t ]] m 
do not include the vertex M ai ^. 

Lemma 12.61 leads to the following decomposition theorem. 
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2. 3. The decomposition process for the polytope S 

Theorem 2.1 (decomposition theorem). If the polytope [[a\, 02, ... , a s ], 

[bi, b 2 , • • • > bt]]m is n °t simplicial, then it can be decomposed into the union 
of two convex polytopes (not always simplicial). The expression is 



Here conv-fS 1 } denotes the convex hull of the set S of points . 
Proof. This follows from Lemma 12.61 



Based on Theorem I2.1[ the polytope S can be easily subdivided into 
simplices. 

Example 1. Show the simplicial subdivision of the following convex poly- 
tope 

L 2 '■= {(Vu - - - , 2/5) I 2/i + 2/2 — 2/3 — 2/4 — 2/5 < 0, (j/i, . . . ,y 5 ) T G A 5 }. 

Solution. Denote L~ as [[1,2], [3,4,%. We know that [[1,2], [3,4,5]] 5 is 
not simplicial by Lemma 12.41 Using Theorem 12.11 we have 



[[1, 2], [3, 4, 5]] 5 = conv{M li3 , [[2], [3, 4, %} U conv{M 1)3 , [[1, 2], [4, %}. 
By Lemma [21] we know that both [[2], [3, 4, 5]] 5 and [[1, 2], [4, 5]] 5 are not 



simplicial. Therefore we repeatedly apply Theorem 12.11 to them and have 



[[2], [3,4,% 

= conv{M 2 , 3 , [[ ], [3,4,%}Uconv{M 2 , 3 , [[2], [4,%}. 

= conv{M 2 , 3 , [[ ], [3, 4, %} [j conv{M 2 , 3 , M 2 , 4 , [[ ], [4, %} \J 

conv{M 2 , 3) M 2)4 ,[[2],[%} 
= conv{M 2j3 , e 3 , e 4 , e 5 } [j conv{M 2j3 , M 2)4) e 4 , e 5 } |J 

conv{M 2i3 ,M 2i 4,M 2i5 ,e 5 }. 

[[1,2], [4,% 

= conv{M 1)4 , [[2], [4,%}Lkonv{M M , [[1,2], [5]] 5 }. 

= conv{M M , M 2j4 , [[ ], [4, %} \J conv{M li4 , M 2j4 , [[2], [%} \J 

conv{M M ,[[l,2],[5]] 5 } 
= conv{Mi j4 , M 2)4 , e 4 , e 5 } [j conv{Mi j4 , M 2)4 , M 2j5 , e 5 } \J 

conv{M M ,M li5 ,M 2! 5,e 5 }. 



[[ai,a 2 , • • • ,a s ), [bx,b 2 , . . . ,b t ]) m 
conv{M ai)6l , [[a 2 , . . . , a s ], [h, b 2 , . . . , b t ]] m } [j 
conv{M aij6l , [[ax, a 2 , . . . , a s ], [b 2 , b t ]] m }. 
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Finally we get the expression of simplicial subdivision of [[1,2], [3,4, 5]] 5 , 
[[1,2],[3,4,5]] 5 

= conv{Mi i3 , M 2 , 3 , e 3 , e 4 , e 5 } \J conv{Mi i3 , M 2j3 , M 2j4 , e 4 , e 5 } |J 
conv{M li3 , M 2)3) M 2i4 , M 2 , 5 , e 5 } U conv{M li3 , M 1)4 , .1/,.,. e 4 , e 5 } U 
conv{M li3 , M 1)4 , M 2i4 , M 2 , 5 , e 5 } U conv{M li3 , M 1>4 , M lj5 , M 2j5 , e 5 }. 

So [[1, 2], [3, 4, 5]]5 is a union of six 4- dimensional simplices. 

We summarize the decomposition process of Example [1] into the following 
algorithm. 

Algorithm 1 (Vmatrix) 

Input: The expression of polytope [[cti, a 2 , . . . , a s ), [bi, b 2 , . . . , b t ]) m . 
Output: Simplices D X) D 2 , . . . , D p (denoted by matrices) such that 

p 

[[ai,a 2 , • • - ,a s ], [&i,6 2 , . . .,b t )) m = (J A, int(A) p)int(D 3 -) = 0, for z ^ j. 

i=i 

VI: Let F := {[[ai, a 2 , . . . , a s ], [6 X , 6 2 , . . . , b t )) m }, temp := 0. 
V2: When F ^ 0, repeat the following procedures 

V21: Choose a polytope iV 6 F . If iV is simplicial, then 

temp := temp U {N}, F :— F\ {N}. 
V22: If the polytope N is not simplicial, then by Theorem [2J] de- 
compose it into two convex polytope Bi, B 2 . 
F := F\ {N} U {Si, B 2 }. Go to step V2. 
V3: Return temp. 



We have written a function in Maple [31] to implement the above algo- 
rithm. 

Lastly, we will present a formula for computing the number of simplices 
given by the polytope subdivision. 

Lemma 2.7. According to algorithm Vmatrix, the convex polytope [[1, ...,&], 
[k + 1, . . . , m]] m (0 < k < m — 1, m > 2) can be subdivided just into f(k, m) 
simplices, where 



f(k,m) = 




{m - 1)! 
k\(m - 1 - fc)!" 
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We know that f(k, m) has the same recurrence formula as binomial coeffi- 
cients by Theorem 12. II Thus the proof of Lemma |2~7T1 is easy via an induction 
argument. This formula will be used to estimate the cost of Algorithm 2 in 
the next section. 



3. Determining algorithm for copositive matrices 



In this section, we will present the complete determining algrorithm of a 
copositive matrix. 

Given annxn symmetric matrix 



A = [ay] 



an oc 
a A 2 



compute A (see (JT])) 



.4 



an ol 
a DA 2 D 



Let B = a\\DA 2 D — aa T , and let 

a = (sign(ai 2 ), • • • , sign(ai„)) T = (&,..., /3„_i) T 

Define the projection operator Proj of the matrix A as follows, 

• If A > 0, i — 1, . . . , n — 1, then 

Proj (A) = {ZM 2 -D}. 

• If there is at least one -1 in fy, then 

Proj (A) = {DA 2 D, W?BW 1} W^BWp}. 

Here the matrices W%, . . . , W^, is fixed by the simplicial subdivision of the 
convex polytope T~ (see (T5])). 

Algorithm 2 (COPOMATRIX) 



Input: Symmetric matrix A £ R nx?1 (r;, > 2). 
Output: A is copositive, or A is not copositive. 
CI: Let F := {A}. 
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C2: Repeat the following steps for the set F. 

C21: If the set F is empty, then return U A is copositive". 

C22: Check the (l,l) th entry of every matrix K in set F . If at least 

one of them is negative, then return Ll A is not copositive". 
C23: Compute the projective set P := \J K f Pto](K) of set F. 

F := P\ {the nonnegative matrices of P}. Go to step C21. 



Note that the above algorithm is also valid for 2x2 matrices. Fur- 
thermore, for strictly copositive matrices we can also formulate a similar 
algorithm. 

The correctness of the algorithm COPOMATRIX is guaranteed by Lemma 
II. H and the algorithm obviously terminates. The cost of the algorithm 
mainly depends on the number of simplicial subdivisions of the polytope. 
According to Lemma 12. 7| we can estimate that in the worst case it is at 
most: 

<(^) +1 )<(^) +1) ...((9 +1 ) 

< (2 n - 3 )(2 n ~ 4 )---(2)(2) 

_ o(n-2)(n-3)/2+l 

The bound 2^ n ~ 2 ^ n ^ 3 ^ 2+1 is already much lower than doubly-exponential 
cost of CAD ' 2,9 l We have written a function in Maple to implement the 
algorithm COPOMATRIX. For non-commercial request, we will offer for 
free. Please sent e-mail to the address 

yaoyong@casit.ac.cn, or, 

j.jia.xu@gmail.com. 
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